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Abstract. We propose a new complexity measure of space for the BSS 
model of computation. We define LOGSPACEvk and PSPACEvk complex- 
ity classes over the reals. We prove that LOGSPACEvk is included in 
NCg n Pw, i-e. is small enough for being relevant. We prove that the 
Real Circuit Decision Problem is Pu-complete under LOGSPACEvk re- 
ductions, i.e. that LOGSPACEm/ is large enough for containing natural 
algorithms. We also prove that PSPACEvk is included in PARr. 
Keywords: BSS model of computation, weak model, algebraic complex- 
ity, space. 



Introduction 

The real number model of computation, introduced in 1989 by Blum, Shub and 
Smale in their seminal paper ,BSS8£,, has proved very successful in providing 
a sound framework for studying the complexity of decision problems dealing 
with real numbers. A large number of complexity classes have been introduced, 
and many natural problems have been proved to be complete for these classes. 
A nice feature of this model is that it extends many concepts of the classical 
complexity theory to the broader setting of real computation; in particular a 
question Pr ^ NPr has arisen, which seems at least as difficult to prove as the 
classical one, and several N PR-complete natural problems have been exhibited. 

It has been soon pretty obvious, however, that all features of the classical 
complexity theory could not be brought to this setting. In particular, the only 
complexity measures considered so far were dealing with time, and not, say, 
space: in 1989, Michaux proved in IMic89| that, under a straightforward notion 
of space, everything is computable in constant space. Therefore, no notion of 
logarithmic or polynomial space complexity exists so far over the reals. A way to 
deal with this situation has been to define parallel complexity classes in terms 
of algebraic circuits, such that the NCr and the PARr classes. 

This model of computation has also long been criticized for being unrealistic: 
the assumption that one could multiply two arbitrary real numbers in constant 
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time was the usual target, that Koiran faced in |Koi97| by defining a notion of 
weak cost that increases the cost for repeatedly multiplying or adding numbers. 

Inspired by his approach, we propose here a new measure of space for the 
real number model, denoted as weak space, such that a repeated sequence of 
multiplications or additions on a number increases its size. Our notion allows 
us to define a logarithmic space complexity class, that falls within NC^. We 
also prove that this class is large enough for containing natural algorithms: in 
particular, we prove a PR-completeness result under LOGSPACEvy reductions. 

The paper is organized as follows: in Section ^ we recall concepts and no- 
tations from the BSS model of computation. We define machines, circuits, and 
some major complexity classes. In Sectional we briefly recall Michaux's result, 
and sketch a proof. In Section |31 we briefly introduce Koiran's notion of weak 
cost, and state some of the major results related to this notion. Then, we intro- 
duce our notion of weak size in Section 01 and state our results. 

1 A Short Introduction on the BSS Model 



In this section, we list the notations used in the paper, and recall some 
basic notions and results on the BSS model. A comprehensive reference for these 
notions is [HnSHMl. 

1.1 Notations 

For an integer c £ Z we define its height as [log(|c| + 1)] . The height of an integer 
is the number of digits of its binary encoding. We also define M* — IJ^gpjM". 

1.2 Real Machines 

We consider BSS machines over R as they are defined in jBSS89IBCS"S98) . 
Roughly speaking, such a machine takes an input from R* , performs a number 
of arithmetic operations and comparisons following a finite list of instructions, 
and halts returning an element in R* (or loops forever). Such a machine can be 
seen as a Turing machine over R. It essentially consists in a finite directed graph, 
whose nodes are instructions, together with an input tape, an output tape, and a 
bi-infinite work tape, equipped with scanning heads. The instructions can be of 
the following types: Start, Input (reads an input value). Output (writes an out- 
put value). Computation (performs one arithmetical operation on two elements 
on the work tape), Constant (writes a constant parameter Ai G R), Branch 
(compares two elements, and branches accordingly). Shift, Copy and Halt. 

For a given machine M , the function ipM associating its output to a given 
input X e R* is called the input- output function. We say that a function / : 
R* — > R* is computable when there is a machine M such that / = fM- 

Also, a set L C R*, or a language is decided by a machine M if its charac- 
teristic function xl:R*^{0,1} coincides with ipM- 



This model of computation allows one to define complexity classes. In partic- 
ular, Pr is the set of subsets of R* that are decided by a real machine that works 
in deterministic polynomial time. Similarly, NPr is the set of subsets of M* that 
are decided by a real machine that works in nondeterministic polynomial time 
i.e, a; G M* is accepted if and only if there exists j/ G R* of polynomial size such 
that the machine accepts {x,y). 

1.3 Configurations 
Definition 1. Configurations 

— A configuration of a machine M is given by an instruction q of M along 
with the position of the heads of the machine and three words Wmput G ^* , 
Wwork e R* • U'ouLpuL "= R* that givc the contents of the input tape, of the 
work tape and of the output tape. 

— A transition of a machine M is a couple {ci,Cj) of configurations such that, 
whenever M is in configuration Ci, M reaches cj in one computation step. 

Definition 2. Configuration Graph 

For a given machine M and a given set C of configurations of M , we define the 
configuration graph of M onC to be the directed graph with vertexes all elements 
in C, and edges all transitions of M between elements in C. 

1.4 Algebraic Circuits 

We introduce the notion of algebraic circuits, that allows to denote parallel 
computations and to define complexity classes below Pr. 

Definition 3. Algebraic Circuit 

An algebraic circuit C is a sequence of gates (Gi, . . . , G^) of one of the following 
types: 

1. Input gates: d = Xi, takes the input Xi from R, 

2. Arithmetic gates: perform the operation * to the outputs of gates Gj and Gi, 
j,l <i and * G {+, -, .,/}, 

3. Constant gates: Gi = Ai, Ai G M, 

4- Sign gates: If Gj > then Gi = 1 else Gi =0, j < i. 

If a circuit has n input gates, we can suppose that they are the first ones, 
Gi , . . . , G„ . // moreover the last node G^ is a sign node, we shall say that C 
is a decision circuit. 

An algebraic circuit is a finite directed graph with no loops: its size is the 
number of gates, and its depth is the length of its longest path, starting from an 

input gate. 

Algebraic circuits extend the classical notion of boolean circuit to the BSS 
setting, and allows one to define the NCr hierarchy of complexity classes: 



Definition 4. NCr 

For all i £ N, NCr is the class of real decision problems decided by a P-uniform 
family of circuits of polynomial size and of depth bounded by C'(log*(n)), and 

NCr = U NCjj. 

As remarked by Poizat in jPoi95| . in the definition above the uniformity of 
the family can be considered relative to the classical Turing model. Hence, a 
P-uniform family of algebraic decision circuit is such that there exists an finite 
enumeration of all the constant gates in the family. There exists then a P time 
Turing machine which, on input n, k, outputs a discrete description of the k^^ 
gate of the 71*'' circuit of the family. 

Proposition 1. \Cuc92}/ 

NCr C Pr. 



2 Michaux's Result 

This section is devoted to a brief exposition of Michaux's Result |Mic89j , which 
states that a straightforward measure of space fails in differentiating one algo- 
rithm from another. In this section, we will use the following notion of space as 
a complexity measure: 

Definition 5. Unit Space 

Let M be a machine over M, and let c be a configuration of M . We define 
USize{c), the unit size of c to be number of non-empty cells on the work tape 
at configuration c. Assume that on an input (xi, . . . , Xn), the computation of M 
ends within t computation steps. The computation follows a path co, . . . , Q. We 
define the unit space used by M on input (xi, . . . , Xn) to be 

U Space{M, (xi, . . . , Xn)) = max U Size{ck)- 

Assume that the running time of M is bounded by a function t. We define the 
unit space used by M on input size n to be 

U Space{M, n) — max U Space{AI, (xi, . . . , Xn))- 

(xi,...,Xn)^R^^ 

This notion of unit space is essentially the same as the classical notion of space 
for Turing machines. While in the classical Turing model this notion gives rise 
to a whole hierarchy of complexity classes like LOGSPACE, PSPACE, interlaced 
with the time hierarchy, this is not the case in the real setting. In order to precise 
a bit how unit space behaves on the reals, let us begin with the following well 
known technical result. 



Lemma 1. Let M be a real machine with parameters Ai, . . . , A,„, whose run- 
ning time is bounded by a function t. Let n G N. On any input xi, . . . G R", 
at any computation step k < t{n), any non-empty cell on the work tape, say 
ei, contains the evaluation of a rational fraction fi^k G "^{Xi, . . . , Xn+m) on 

(xi , . . . , Xn , ^1 1 ■ ■ ■ 7 ) • 

Proof. Details arguments can be found in |Mic89IPoi95IKoi97| . We only sketch a 
proof here. The key argument is that, at any computation step k, the content of 
any cell e; is obtained from the input values and the parameter values by a finite 
sequence of arithmetical operations. Therefore, the value in e; is the evaluation 
of a rational fraction fi^k £ . . . , X„+rn) on (xi, . . . , x„, ^i, . . . ,Am)- 

Proposition 2. ]Mic89f Let L C K* by a real language decided by a machine 
M in time bounded by a function t. There exists a constant fc G N and a machine 
M' deciding L in unit space k. 

Proof. We only sketch the proof here. The interested reader can find more ex- 
planations in |Mic89IPoi95] . 

Rational fractions with integer coefficients can be easily encoded in binary, 
therefore, by Lemma^ any configuration of M can also be encoded in binary. It 
suffices to realize that this binary encoding can be embedded into the digits of 
only two real numbers. Then, there exists a machine M' simulating M with only 
a constant number of real registers, among which two are needed for encoding 
the configurations of M. 

3 The Weak BSS Model by Koiran 

3.1 Definitions 
Definition 6. Weak Cost 

Let M be a machine whose running time is bounded by a function t, and let 
Ai, . . . , Am be its real parameters. On any input xi, . . . ,Xn, the computation of 
M consists in a sequence cq, . . . , ct,t < t{n) of configurations. To a transition 
Ck,Ck+i in this sequence we associate its weak cost as follows: 

— If the current instruction of c^ is a computation node, let ci be the current 
cell on the work tape: the transition Ck, Ck+i consists in the computation of a 
rational fraction fi+i^k+i = 5i+i,/c+i/^i+i,fc+i G Z( 

which is placed on the cell e;+i in Ck+i. The weak cost of the transition 
Ck,Ck+i is defined to be the maximum of deg(gi+i^k+i), deg(hi+i^k+i), and 
the maximum height of the coefficients of gi+i^k+i and hi+i.k+i. 

— Otherwise, the weak cost of the transition Ck,Ck+i is defined to be 1. 

The weak running time of M on input xi, . . . ,x„ is the sum of the weak costs 
of the transitions in the sequence Cq, . . . , Cj. 

The weak running time of Af is the function that associates with every n the 
maximum over all x G M" of the running time of M on x. 



3.2 Some Results 



Lemma 2. [KoiQUl ^ function is polynomial-time in the weak BSS model if 
and only if it is polynomial-time computable in the standard BSS model and 
the rational fractions fi k have polynomial degree and coefficients of polynomial 
bit-size. 

Let Pw (respectively NP^f) be the set of real languages decided in deter- 
ministic (resp. nondeterministic) weak polynomial time, and EXP^/ be the set 
of real languages decided in weak exponential time by a real machine. 

Proposition 3. 

NCr Ph/ C Pr C NPw = NPr C PARr C EXPw 

Pw- C NPvi/ = NPr is from |(JSS94| . NPr C EXPw from |Koi97| (where it is 
shown that the inclusion is strict). The missing items can be found in [BCSS98J. 

4 Weak Size and Space 
4.1 Definitions 

Instead of considering a unit size for all values on the work tape, which allows 
one to decide every decidable language in constant space, we would like to have 
a notion of size for the values computed on the work tape. The weak size of a 
computed value is a reasonable upper bound for the size of a boolean description 
of the corresponding rational fraction with integer coefhcients. The weak size of 
a configuration is then the sum of the weak sizes of all computed values on the 
work tape in this configuration. 

Yet, we need to precise a bit more the idea. Our purpose is to have a "nice" 
measure of space, allowing one to define a reasonable logarithmic space class. 
A trivial rational fraction like fi{Xi, . . . , X„, Ai, . . . , Am) — Xi has clearly a 
boolean description of size 1, while, for describing fn{Xi, . . . , X„, Ai, . . . , A„i) = 
Xn, one would need [log(ri + l)] digits (for encoding the variable index). It seems 
rather unsatisfactory that a logarithmic space configuration may have a loga- 
rithmic number of occurrences of /i, but only a constant ones of /„. This feature 
can be corrected by allowing a permutation of the input variables, provided the 
permutation is simple enough, i.e can be described in logarithmic boolean space. 
In this paper, we have restricted ourselves to circular permutations, that can be 
described by an offset in {0, . . . ,n — 1}. 

Definition 7. Weak Size 

Assume Ai, . . . , Am G are given real numbers, and let g e . . . , Xn+m] 

be a real polynomial with integer coefficients. Define a real polynomial gAi,...,A^ — 
g[Xi, . . . , Xn, Ai, . . . , Am], with free variables Xi, . . . , X„. Let 0<O<n, OsN 
be a number, the offset. To g, Ai, . . . , Am £ and 0, we associate the follow- 
ing: 



— deg{g) is the degree of g. We will write D{g) for \log{deg{g) + 1)]. 

— VarAi....^A^{g) C {Xi, . . . , Xn} is the set of input variables on which 5^1,. 

effectively depends. 

— Rai....,a„,,o{9) = max{z + mod n} for Xi G VarAi,...,A^{g), is the range 
of g. We will write R{g) for |"log(i?Ai,..,A„,o(5) + 1)1 • 

— N{g) £ N is the number of non-zero monomials of g. 

— S{g) € Z is the maximal absolute value of the integer coefficients of g. We 
will write S{g) for \\og{2S{g) + 1)]. 

— VcAi,...,Am{9) is the maximum, for every monomial of g, of the number 
of input variables on which it effectively depends. We will write y{g) for 
VcA^,...,Ar^{g)- 

The weak size SA-i_,...,Am,o{g) of g is defined as follows: 



For a rational fraction f = g/h we take the weak size of f to be the maximum 
of the weak sizes of g and h. 

It is clear that the weak size of g bounds the size of a boolean encoding of g, 
where g is presented as a sum of monomials modulo a circular permutation of 
the variable indexes. We do not take into account succinct boolean descriptions 
of factorized polynomials to ensure the tractability of our measure. 

This measure of size for an element on the work tape naturally yields a notion 
of weak space for the given work tape, as follows: 

Definition 8. Weak Space 

Let M be a machine with real parameters Ai, . . . , Am- Let Ck be a configuration 
of M, with the corresponding input xi,. . . ,Xn G M" . We define: 

— ei,...,ej to be the non-empty part of the work tape in the configuration Ck- 

— For any non-empty cell e/ in Ck, we denote by fi^k G 'Z{xi, ■ . ■ , Xn, Ai, . . . , A^) 

the rational fraction it contains. 

The weak size of the work tape at the configuration Ck is then: 



Assume that the running time of M is bounded by a function t. For a given 
input xi, . . . , Xn, the computation of M consists in a sequence cq, ■ . ■ ,Ct,t < t{n) 
of configurations. 

The weak running space of M on input xi,. . . ,Xn is the maximum for all 

configurations ca, . . . ,ct of their weak size. 

The weak running space of M is the function that associates with every n 
the maximum over all x of the running space of M on x. 



SA,,...,A^,o{a) = N(5) (S(3) + \/{g).R{g) + V(c;).D(c;)) 



(1) 




Definition 9. Complexity Classes 

— A language L C R* is in LOGSPACEvf if and only if there exist a machine 
M and a constant fc G N such that, for all n N, on input x £ R", M 
decides whether x L in weak space less than fclog(n). 

— A language L CM.* is in PSPACEvi^ if and only if there exist a machine M 
and two constants k,d G N such that, for all n E N, on input x G M", M 
decides whether x E L in weak space less than kn'^ . 

— A function / : R* R* is in FLOGS PAG if and only if there exist a 
machine M and two constant k,m £ N such that, for all n G N, on input 
X G R", and computation cq, . . . ,Ct of M on x: 

1. M computes f{x) in weak space less than fclog(n). 

2. for every configuration Ci with current node an output node and current 
cell ei, the weak size of the content of ei is less than m. 

In the definition of FLOGSPAGEvf, the output consists in a sequence of real 
values of constant weak size in the input. This ensures that one can compose 
FLOGSPAGEvF algorithms, and that the result of the composition remains an 
FLOGSPAGEvF algorithm. This is necessary for defining notions like logarithmic 
space reductions and for obtaining completeness results. 

4.2 What Michaux's Result Becomes 
Lemma 3. There exists L C R* such that: 
~ LgPw 

— for all k E N, L is not decidable in weak space less than k. 

Proof. Let p{Xi, . . . , X„) — Xi + . . . + Xn, and consider the set L of points 
(xi, . . . , Xn) G R", such that p{xi, . . . , Xn) equals 0. Assume L is decided by a 
machine M . It is well known that the set of inputs accepted by a BSS machine 
is semi-algebraic, therefore, L can be described as a finite union of sets given by 
systems of polynomials inequalities of the form 

s t 

/\ F,{Xi, . . . , X„) = A /\ Gj(Xi, . . . , X,0 > 0, 
1=0 i=o 

where the values Fi(xi, . . . ,Xn) and Gj(xi, . . . , Xn) are effectively computed by 
M . Since L has dimension n — 1, at least one of these sets must have dimension 
n — 1. Since the set described by the G'jS is open, it must be nonempty, and 
then it defines an open subset of R". All the polynomials Fi vanish on that 
nonempty open subset of L. Since this open subset of L is clearly infinite, and 
p is an irreducible polynomial, all the polynomials Fi must vanish on the whole 
set L. It is then a well known result f BCSS98 . Proposition 2 p. 362) that the 
polynomials Fi are multiples of p. Also, at least one of these Fi is a non-trivial 
multiple of p. 

It is clear that p{xi, . . . ,Xn) has weak size at least nlog(n), and so does this 
non-trivial multiple of p. Therefore, M decides L in weak space at least n log(n) . 



4.3 Structural Complexity Results 



Theorem 1. 

LOGSPACEm/ CP^n NCr, 
PSPACEvK C PARr. 

Proof. In a first step, we prove LOGSPACEvf C Pr. The key argument is an 
upper bound for the number of configurations of weak size at most fclog(n). 
Consider a machine M, with t nodes. For a fixed input size n, and an offset 0, 
simple counting arguments show that the number of rational fractions of weak 
size at most B for some i? e N is bounded by , for some a g N. It follows 
that the number of possible work tape contents of weak size B, for the same 
fixed offset, is bounded by {2a'^)^ . Taking into account all possible values for 
the offset, the scanning head positions and the current node of the machine, the 
number of configurations of weak size at most B is then bounded by tn?B{2a^)^ . 
When B = k\og{n), this bound is polynomial. 

LOGSPACEvi/ C Pw follows then by LemmaEl since all rational fractions of 
logarithmic weak size have clearly polynomial degrees and coefficient heights. 

LOGSPACEh' C NCr is then proven along the lines of |Bor77j : given a 
LOGSPACEvF machine M, we exhibit a NCr construction of its configuration 
graph. This construction involves some numeric computation, in order to check 
whether two given configurations are connected, and produces a boolean de- 
scription of the configuration graph of Ad . Next, it suffices to decide whether 
the input and accepting configurations are connected in this graph: this is the 
classical reachability problem, which is decidable in the boolean class NC^. 

PSPACEw C PARr is a corollary 

4.4 Completeness Results 

Definition 10. \CT9^ Real Circuit Decision Problem (CDPr) 

Input: (C, x), where C is an arithmetic circuit with k input gates and a; G M*^. 

Question: Does C output 1 on input x? 

It has been shown in |CT92| that CDPr is Pg-complete under NCR-reductions. 

Theorem 2. CDPr is Pfi-complete under VLOGSP IKC^w -reductions. 

Proof. The proof follows |CT92j . The reduction happens to be in FLOGSPACEvf- 

We have stated this completeness results under FLOGSPACEvi/ reductions. By 
Theorem ^ it is clear that FLOGSPACEvk reductions are in Pw n NCr. The 
problem considered has already been proven complete under first-order reduc- 
tions |CM96j . which also happen to be in Pw n NCr. Yet, it remains unclear 
how the two types of reductions compare. 



5 Concluding Remarks and Open Questions 



In the discrete model, space has proven to be a very relevant complexity measure. 
Many natural problems have been found in LOGSPACE, and many others in NC^ 
whose membership in LOGSPACE is unclear. Wc believe that weak space may 
play the same role in the real setting. An argument in this direction is the 
following remark: consider a real algorithm that reads an input, normalizes it to 
{0, 1} with some step function, and applies a boolean LOGSPACE procedure. Real 
complexity analysis until now only allowed one to say that such a real algorithm 
belongs to Pw H NCg: the algorithmic flavor behind it was lost. However, it is 
now clear that such an algorithm belongs to LOGSPACEvf- An important task 
now is to exhibit some natural problems in LOGSPACEvk- Others in NC^ or Pw, 
not easily in LOGSPACE^y, may also be of interest. 

Structural results remain also to be found. In particular, it needs to be 
checked whether the following conjecture holds: 

Conjecture 1. 

NC^ L0GSPACEv^^, 
LOGSPACEvK C NC4 ^ LOGSPACE C NC^ 
Similar questions arise also for PSPACEv^'. 
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